Calculate_Gauss_Seidel_preconditioner Function

public function Calculate_Gauss_Seidel_preconditioner(A) result(L)

Arguments

Type IntentOptional Attributes Name
real(kind=dp), intent(in), DIMENSION(:, :) :: A

Return Value real(kind=dp), DIMENSION(SIZE(A, 1), SIZE(A, 2))


Calls

proc~~calculate_gauss_seidel_preconditioner~~CallsGraph proc~calculate_gauss_seidel_preconditioner Calculate_Gauss_Seidel_preconditioner proc~diag Diag proc~calculate_gauss_seidel_preconditioner->proc~diag

Source Code

    FUNCTION Calculate_Gauss_Seidel_preconditioner(A) RESULT(L)
        REAL(dp), DIMENSION(:, :), INTENT(IN) :: A
        REAL(dp), DIMENSION(SIZE(A, 1), SIZE(A, 2)) :: L
        INTEGER :: N, i, j

        N = SIZE(A, 1)

        L = 0.d0

        IF (ANY(Diag(A) < epsi)) STOP "ERROR :: Zero diagonal in Gauss-Seidel preconditioner"
        FORALL(i=1:SIZE(A,1), j=1:SIZE(A,2), i >= j) L(i, j) = A(i, j)

    END FUNCTION Calculate_Gauss_Seidel_preconditioner